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The inhomogeneous 3-Kelvin (3K) phase of the eutectic Sr2Ru04 with Ru inclusions nucle- 
ates superconductivity at the interface between Ru and Sr2Ru04. The structure of the interface 
state and its physical properties are examined here. Two superconducting phases are identified 
between the transitions to the bulk phase at 1.5 K and to the 3K phase. The nucleation of 
the 3K phase results in a state conserving time reversal symmetry, which generates an intrin- 
sically frustrated superconducting network in samples with many Ru inclusions. At a lower 
temperature (> 1.5 K), a discontinuous (first order) transition to an interface state breaking 
time reversal symmetry is found leading to an unfrustrated network phase. It is shown that 
this phase transition located at a temperature between 1.5 and 3 K would yield the anomalous 
property showing that the critical current in such a network depends on the sign of the current, 
reproducing recent experimental observations. 
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1. Introduction 

Sr2Ru04 is a quasi-two-dimensional strongly corre- 
lated metal showing unconventional superconductivity 
with a bulk critical temperature — 1.5 K.^'^ There 
is strong evidence that the superconducting phase has 
a spin-triplet odd-parity character and breaks time re- 
versal symmetry. The most likely order parameter has a 
chiral p-wave symmetry corresponding to a pairing state 
with the orbital structure ± iky and the triplet spin 
configuration corresponding io Sz — 0. This state is 
usually represented by the vector gap function d{k) = 
z{kx ± iky) and is described as a two-component order 
parameter, for example as d{k) — rj^zkx + rjyzky. 

Some years ago, an intriguing feature of Sr2Ru04-Ru 
eutectic samples was discovered. Some of the Ru seg- 
regates into ^m-sized inclusions embedded within the 
parent material Sr2Ru04. These samples show an onset 
to inhomogeneous superconductivity at approximately 
T* « 3 K, which turns into the bulk superconduct- 
ing phase at T = T^^.^''^ This higher-temperature phase 
(-rpb < T < T*) has been called the "3-Kelvin" phase 
(3K phase). The early suggestion that the superconduc- 
tivity in the 3K phase has filamentary nature, nucleating 
at the interface between Ru-inclusion and Sr2Ru04,^ re- 
ceives strong support on the basis of the behavior of the 
upper critical fields Hc2, which is enhanced and shows a 
characteristic sublinear dependence on \T — This 
theory is based on the assumption that, for some as- 
yet unknown reason, the pairing interaction is enhanced 
in Sr2Ru04 in the vicinity of the interface.^ Interest- 
ingly, the effective symmetry lowering at the interface 
implies the nucleation of a superconducting state, which 
is different from that of the bulk state of Sr2Ru04 and 
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does not violate time reversal symmetry. For the two- 
dimensional order parameter space of the p-wave state, 
i.e., d(k) — rjxzkx -\-r]yzky, quite general arguments lead 
to the conclusion that the pairing state with the orbital 
symmetry k ■ n — (the p-wave lobe parallel to the in- 
terface) is realized at T = T* (n is the normal vector of 
the interface).^ Consequently, the transition from the 3K 
phase to the bulk superconducting (chiral p-wave) phase 
is not merely a percolation transition as in the case of 
conventional inhomogeneous superconductors. At least, 
time reversal symmetry breaking (TRSB) has to occur 
on the way to the bulk phase. 

In this study, we investigate the evolution of the 
3K phase towards the bulk superconducting phase of 
Sr2Ru04 in the temperature range < T < T* . In this 
context, we will discuss three different phases for a closed 
interface of a Ru inclusion, which we call the A-, A'- 
and B-phases. The A-phase appears at T* and is time- 
reversal-conserving. At lower temperatures, a transition 
to a time-reversal-symmetry- violating phase named A'- 
phase (B-phase) occurs, which may have the same (dif- 
ferent) topology as the A-phase. We will show that the 
A- and A'-phases introduce phase frustration, if we con- 
sider a network formed by the superconducting interfaces 
states of many Ru inclusions.^ On the other hand, the 
topology of the B-phase yields a nonfrustrated network. 
Our investigation shows that the additional transition is 
of first order from the A-phase to the B-phase; it breaks 
time reversal symmetry and simultaneously changes the 
topology of the states on an inclusion. The B-phase is 
then topologically identical to the bulk phase, which is 
eventually reached by percolation. The additional transi- 
tion is accompanied by observable effects such as charac- 
teristic features of quasiparticle tunneling^"^*^ and critical 
current. 
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2. Model Formulation 

To illustrate the most relevant features of the 3K 
phase, we first consider a single Ru inclusion modeled as 
a cylinder of radius R whose central axis lies along the 
z axis of the tetragonal crystal lattice of Sr2Ru04. Here, 
we focus on the superconductivity in the two-dimensional 
x-y plane and use the in-plane p-wave spin-triplet pairing 
as the dominant superconducting instability. Assuming 
for the bulk phase the chiral p-wave spin-triplet state, we 
represent the order parameter by the vector gap function 
as d{r, k) = rj+{r)z{kx + iky) + r]^{r)z{kx — iky). With 
this order parameter, we now write a Ginzburg-Landau 
(GL) free energy for the region belonging to Sr2Ru04 
around a single cylindrical Ru inclusion (r > i?), 

F = foIr>Rd'r [r{r)\r,\^ + l{\r,\^ + 2\rj+\'\n.\^} 



+ {|£>r/+|2 + |£)r?_|2} + 2k2(V x A)^ 
m{D_T,+)*{D+n_)+c.c.}]. 



(1) 



We use dimensionless units for our formulation, and the 
notation t] = {r]+,r]-), T{r) = T/Tc{r) - 1. With the 
r-dependence of the critical temperature T'c(r), we in- 
corporate the local enhancement of the pairing inter- 
action. We introduce the operators D = V — iA with 
D± = Dx±iDy and A as the dimensionless vector poten- 
tial. The units are chosen so that the zcro-tcmpcraturc 
coherence length is the unit length, the order parame- 
ter reaches or |t?-|^ = — r in the uniform bulk state 
below (2-fold degenerate state), and the vector poten- 
tial is given in units of <l>o/27r^o (^o = hc/2e is the flux 
quantum). Moreover, k is the usual GL parameter and 
/o is the free energy per \mit volume. A straightforward 
variational calculation shows that the coordinate separa- 
tion of the order parameter in eq. (1) is adequate here 
with N as the qiiantiim mimbcr characterizing different 
topological sectors of the order parameter (see eq. (3)). 

In order to model the narrow region of enhanced super- 
conductivity at the interface, we introduce the spatially 
dependent Tc(r) in the following form: 

Tc(r) = tI; + To/ cosh[(r - R) /d] for r > R. (2) 

Here, d is the width of the region close to the interface 
where Tc is locally enhanced with a maximum Tq at the 
interface. At this moment, we ignore the magnetic field 
{A = 0) and investigate first the nucleation of supercon- 
ductivity. 

With the given geometry and order parameter choice, 
it is advantageous to turn to cylindrical coordinates 
{r,6,z). The GL equations allow to separate the depen- 
dences of the order parameters in radial and azimuthal 
coordinates into the form (assuming homogeneity along 
the 2-axis), 

»7±(r)=^7±(r)e^(^^l)^ (3) 

where N is an integer corresponding to the phase wind- 
ing number of d{r, k) around the cylinder, since the 
the order parameter is single- valued. This representation 
reflects the cylindrical symmetry of our model geome- 
try and gives rise to the following symmetry properties. 
One finds d dexp{iN'y) for a rotation of the sys- 
tem, 6 ^ 6 + •y and Ok ^ Ok + 7, owing to the factor 



kx ± iky = \k\ exTp{±i0k) and the phase factor of cq. (3). 
Here, we always keep the orientation of the d vector fixed 
parallel to the z-axis. The formulation of the boundary 
conditions is most conveniently performed with the ra- 
dial and azimuthal components of the order parameter 
defined as 



1 



Vr{r,0) = ^{r^+{r)+rj_{r))e^''<' 
ry,(r,^) = -^(7?+(r)-7?_(r))e*^« 



(4) 



These two order parameters are differently affected by 
boundary effects due to their interference effects under 
surface scattering. These effects can be implemented by 
using the following standard boundary condition at the 
interface {r = R); 



dr 



(5) 



where n stands for "r" or "0" . Here, £^ is the so-called ex- 
trapolation length characterizing the boundary effect for 
each component (for example, see, ref. 12). We assume 
— > oo, since the tangential component of a Cooper pair 
wave function at an interface is not affected by (specu- 
lar) surface scatterings.^^ On the other hand, there is 
a suppression of the radial component (sign change or- 
der parameter under reflection at the interface) requiring 
a finite value of £r- In our model, we neglect the su- 
perconducting component induced in the Ru inclusion, 
assuming an interface of very low transmissivity, con- 
sistent with our assumption of the boundary condition. 
This model will now be used to derive a phase diagram 
of the 3K phase. 

3. Phases and Phase Diagram 

For given N and T, we can now determine the r de- 
pendence of the order parameter by varying the GL free 
energy including the boundary conditions. We find three 
states which are of interest for the discussion of the 3K 
phase, which we named the A-, A'-, and B-phases. They 
are distinguished by time reversal symmetry and the 
phase winding number N (sec Fig. 1). The A-phasc is 
a state conserving the time reversal symmetry {rje ^ 
and rjr = 0) and with N = 0. Both the A'- and B- 
phases violate the time reversal symmetry with A'^ = 
and A'' = ±1, respectively {rje ^ and j}r 7^ 0). 

At T = T*, the A-phase nucleates, corresponding to 
?7+(r) = —rj^{r) oc rio{r) according to eq. (4). As a result 
of the boundary conditions {£$ = 00), the order parame- 
ter r]g (r) is the largest at the interface r = R and falls off 

I 1—1/2 

exponentially on the length scale ^'(T) ^ |T(r ^ R) | 
with increasing r. Because N = 0, the internal phase 
structure of the p-wavc pair wave function yields the 
phase difference tt across a Ru inclusion for the order 
parameters ?7±(t) in cq. (3) (sec Fig. 1). 

The A'-phase is also in the topological sector of = 0. 
A second-order transition from the A-phase may lead to 
this phase at a temperature below T*. This transition 
breaks the time reversal symmetry introduced by the 
continuous appearance of the subdominant component 
r]r. From eq. (4), we find that the A'-phase breaks the 
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Fig. 1. Topological structure of the A-, B-, and metastable A'- 
phases. (a) The A-phase corresponds to a time reversal symmet- 
ric state with = (winding property). The tangential com- 
ponent of a pair wave function is arranged in such a way as to 
keep the phase constant around the interface. However, it yields, 
a diametral phase shift tt due to the internal structure of the 
tangential p-wave state, as indicated in the figure. The perpen- 
dicular component vanishes, (b) The A'-phase with Af = is 
TRSB due to a finite perpendicular p-wave component, (c) The 
B-phase is TRSB with A'^ = ±1 and has a topological structure 
compatible with that in the bulk phase, z{kx it iky). 

balance between the amplitudes |?7-i-(»')| and |?7-(r')| be- 
cause ?7r 7^ and rje 7^ 0. In this way, one of the two order 
parameters ri±{r) becomes predominant as the temper- 
ature is lowered, and will eventually decide the chirality 
of the bulk superconducting phase. 

In contrast to the A and A'-phases, the B-phase is a 
state with iV ^ 0. We distinguish the two degenerate 
states 7^ and — with = 1 and 77+ = and 
r]_ ^ with A^ = — 1. Consequently, the dependence on 
cancels in the representation of the order parameter 
in Eq. (3), unlike the A- and A'-phases. Therefore, this 
state has the same symmetry and phase topology as the 
homogeneous bulk phase of Sr2Ru04, as can be easily 



I 1 



-0.: 





FRSB - 








■TRSHf 














\ 


B 




R 




A 


% 












1s1 


t-order 


d=0.5 






tra 




>ILIUII 


R=10— ^ 
R=2 -- X 










f. IK X ill -»-^ 


1 1.5(=T^) : 




[' 3 



Fig. 2. Temperature dependence of the free energy difference 
Fj^—q{T) — Fpf—i{T) between the states with winding numbers 
N = and A = 1. The two plots correspond to those in the 
cases of different-sized Ru inclusions with radii R = 2 and 10 in 
4o units. The first-order transition from the A-phase (N = 0) to 
the B-phase (A = 1) occurs at T' « 2.4 K. 

deduced from Fig. 1. Therefore, the B-phase naturally 
connects to the bulk phase occurring at T^!. 

In the next step, we examine the sequence of tran- 
sitions in the temperature range between T* and Tj^ 
(Fig. 2). It can be anticipated that below the onset of 
superconductivity at T = T*, a second transition at a 
temperature T = T' leads to either the A'-phase or the 
B-phase. The latter transition (A^B) would be of the 
1st order, since it involves a discontinuous change in the 
winding number [N = Q ^ N — ±1). Iw order to de- 
termine which of the two states is reached, we have to 
compare their free energies. The parameters are set to 
T(? — 1.5 K of the pure Sr2Ru04 and d ~ 0.5 in units 
in eq. (2). Tq is adjusted so as to obtain the nucleation 
temperature for the A-phase to T* 3 K. Note that, on 
a qualitative level, the following results are insensitive to 
the choice of these parameters. 

Our numerical evaluation of the variational equations 
with the given parameters yields a transition from the 
A-phase to the B-phase at T' « 2.4 K (see Fig. 2), which 
is of the 1st order for the reason discussed above. The 
continuous transition between the A- and A'-phases has 
a lower critical temperature (T" « 2.2 K ) and is conse- 
quently not realized. We have confirmed that the transi- 
tion temperature T' shows little sensitivity on the radius 
R of the Ru inclusion as long as i? 1 (length unit ^o)- 
Therefore, in the case of the single Ru inclusion, we con- 
clude that the first-order transition (A— >B) occurs at a 
temperature below the 3K phase nucleation at T = T*. 
The B-phase then evolves into the bulk phase at T = 
extending throughout the whole sample. 

4. Nature of the Two Phases 

In this section, we consider two characteristic proper- 
ties of the A- and B-phases that can help identify the 
two phases experimentally. 

4-.1 Spontaneous currents 

We may physically distinguish the A-phase from the 
B-phase by the fact that the latter carries a spontaneous 
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(7V + l)<(r) 
V-{r)- 



1)772 (r) 
drj^{r) 



dr dr 
Obviously, this current vanishes in the A-phase where 
iV = and 77+ = — In the B-phase, on the other 
hand, current flows near the interface. Note that the ex- 
pression of the current in eq. (6) and the numerical result 
shown in Fig. 3 does not include the screening effects (i.e., 
the vector potential is not determined self-consistently) . 
The Meissner-Ochsenfeld screening is of minor impor- 
tance here since the the bulk of the material is not super- 
conducting. The currents are rather small and sparse so 
that it would be difficult to directly observe the magnetic 
field generated by those currents, for example, by scan- 
ning SQUID or Hall probes. Moreover, the observation 
of the magnetic fields by means of /xSR zero-field relax- 
ation rate measurements is likely difficult owing to the 
small volume fraction. Below, we will discuss that critical 
current measurements may be viewed as evidence of the 
realization of phase transition below T*, which leads to 
a TRSB state. 

4-2 Coupling between inclusions 

Some eutectic Sr2Ru04-Ru samples contain regions 
with a rather high density of Ru inclusions. In these sys- 
tems, superconducting condensates nucleating at the in- 
terfaces of neighboring inclusions can overlap, if they are 
separated by a length on the order ^ ~ ^0/ \/W\ only. A 
composite of many inclusions can form a superconduct- 
ing network in the 3K phase whose properties are influ- 
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Fig. 3. Spontaneous current around a Ru inclusion in the B-phase 
originating from the TRSB at T = 1.9 K above T<?. The radius of 
the Ru inclusion is i? = 10 and the interface is located at r = 10 
in 4o units. Meissner-Ochsenfeld screening is not included in this 
result. 



current that flows along the interface. This is a conse- 
quence of the TRSB, while the change in the topology 
of the states {N = 0— )-7V = ±l)is less important in 
this context.^ The expression of the supercurrent is eas- 
ily derived from the GL free energy [eq. (1)] by using the 
derivative with respect to the azimuthal vector potential. 
In the present cylindrical geometry, the radial component 
of the supercurrent vanishes for symmetry reasons. The 
azimuthal component is given as 

2 r 



Fig. 4. Inhomogeneous systems with multiple Ru inclusions. 
Frustration occurs in the A-phase (left panel), while no frustra- 
tion occurs in the B-phase (right panel). Spontaneous currents 
in the B-phase have the same circular orientation. 



enced by the structure of the order parameter on the in- 
terfaces. The region of overlapping condensates on neigh- 
boring inclusions can be viewed as a weak link or a super- 
conductor/normal metal/ superconductor (SNS) Joseph- 
son junction. The energetically favored state is realized 
when the overlapping order parameters have identical 
phases on both interfaces, as in a standard Josephson 
junction. For the A-phase condensate, the requirement 
of identical phases leads to the configurations shown in 
Fig. 4 (left panel): the p-wave lobes at the closest point 
align with the same phase (" O-phase" ) corresponding to 
parallel arrows for the two upper inclusions in Fig. 4. 
Such an arrangement leads to the opposite orientation 

, of the arrows winding around these two inclusions.^ This 
isj a consequence of the phase structure of the p-wave 
Cooper pairs. This type of coupling prefers the to — " 

((5gHnding) configuration of the phases of the supercon- 
ducting order parameters and, consequently, gives rise 
to frustration as depicted for the "7r-phase" coupling at 
the lower inclusion in Fig. 4 (left panel) , analogous to the 
case of a triangle of antiferromagnetically coupled Ising 
spins. It has been earlier speculated that this frustration 
would yield a non trivial spatial dependence of the order 
parameter phase and would be visible in the magnetic re- 
sponse of the system, e.g. in ac susceptibility. Thus far, 
such frustration effects have not been clear observed in 
experiments. 

In contrast to the A-phase, the B-phase has the prop- 
erty in which neighboring inclusions with overlapping 
condensates would energetically prefer a configuration 
where both interfaces carry a state of the same chiral- 
ity, as shown in Fig. 4 (right panel). This ferro-like cou- 
pling does not lead to frustration. The spontaneous su- 
percurrents around the inclusions have the same circu- 
lar orientation, corresponding to a "ferromagnetic" cou- 
pling of the orbital moments as the stable configuration. 
Therefore, the transition between the A- and B-phases 
also influences the network properties, removing order 
parameter phase frustrations. In particular, we expect 
that the circular currents could introduce magnetic flux 
in the voids of the network, which may be interesting for 
the response of the system as well as the critical current. 
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5. Signature of the T'-Transition in Critical Cur- 
rent Measurements 

In this section, we address the question of how the 
transition between the A and B-phases can be observed 
experimentally. Since both phases are filamentary, ther- 
modynamic bulk properties such as specific heat would 
unlikely provide a sufficiently large signal that is measur- 
able. Interestingly, the supcrcurrcnts carried by the net- 
work of filamentary condensates turn out to give a probe 
for a transition within the 3K phase. Recent experiments 
actually gave evidence of a qualitative change in the crit- 
ical current at T = 2.3K}^ To explain this finding, we 
assume that the superconducting network percolates and 
can carry a very small but finite supcrcurrcnt through the 
sample. Some essential properties of the current flow in 
such a network are captured by a simple model configu- 
ration composed of one loop formed by superconducting 
paths, as shown in Fig. 5. The important feature of such 
a SQUID-like structure [Fig. 5 (lower panel)] is the mul- 
tiple connectivity that makes the current flow susceptible 
to magnetic fluxes threading the network. 

In our "network" the supercurrent splits into two 
branches, i.e., 1 and 2, corresponding to upper and lower 
rows of Ru inclusions, respectively, in Fig. 5 (upper 
panel). For our SQUID-model, we assume each of these 
two branches as single Joscphson contact with its speciflc 
sin-like current phase relation: 

I^h+l2 = Id sin (j) + Ic2 sin((/- + a) , (7) 

where we characterize the branches by their intrinsic crit- 
ical currents Id and 7c2- These describe the effective cur- 
rent phase relation in the two chains of the Ru inchisions 
in Fig. 5 (upper panel). The phase coherence of the super- 
conducting order parameter is described by the phases (j) 
and a where the latter takes into account the phase dif- 
ference along the two branches due to the magnetic flux 
$ enclosed in the loop, given in a simple effective form 
as 

a = |^($ + L(7i-/2)). (8) 

Note that $ represents the magnetic flux induced not 

only by external magnetic fields, but also by sponta- 
neous currents running along the Sr2Ru04-Ru interfaces 
in the B-phase. These spontaneous currents are not in- 
cluded in the Joscphson currents, as they need not to 
pass between Ru inclusions. The second term represents 
the contribution to the flux due to the Joscphson cur- 
rents running through the network, specifically through 
the two arms of the SQUID structure, with L as an ef- 
fective self-inductance. For simplicity, we assume the two 
arms to be symmetric in geometry such that their contri- 
butions to the induction are identical with opposite signs, 
while the Josephson coupling strengths are different for 
the two arms (Jd ^ Ic2) in order to introduce asymmetry 
necessary to avoid the cancelation of the current-induced 
contribution to a. 

Now we consider the measurement of the critical cur- 
rent in this network, which corresponds to the maxi- 
mal supercurrent that can be transfered from one end to 




Fig. 5. Simple SQUID network model: Upper panel: The super- 
conducting network formed by the condensates around Ru inclu- 
sions, as depicted here, can be viewed as a multiply connected 
system where the supercurrent can branch and flow in different 
paths. Lower panel: The configuration of the upper panel can be 
modeled as a SQUID-like system with two branches, which in- 
corporates the essential features of a network. Note that, in the 
B-phase, spontaneous currents flow around eaxii inclusion with 
the same circular orientation such that a net flux occurs through 
the loop. 

the other. In the time-reversal-invariant A-phase, sponta- 
neous currents and flux $ are absent (<l> = 0). Obviously, 
the maximal current which can flow through the SQUID 
network is independent of current direction, as can be 
easily verified from eqs. (7) and (8), because the reversal 
of current 7 ^ — 7 is implemented by — >■ — ^ leading to 
I^ ^ -I, {i = 1,2). 

In contrast, in the B-phase, spontaneous supercurrents 
at the interfaces of Ru inclusions generate a finite flux 
$ in the loop, which remains basically unaffected by the 
small Josephson currents. Thus, the analysis of our model 
shows an explicit symmetry breaking for the current re- 
versal 7 —7. The change in the sign of docs not lead 
anymore to a reversed current in the SQUID, because 
the flux $ is flxed in the TRSB phase. Together with the 
contribution of the current induced flux, this leads to 
the effect that the maximal current 7 becomes direction- 
dependent. $ determines the magnitude of the critical 
current difference that changes sign at $ = n$o/2, as 
shown in Fig. 6. 

In their experiment, Hooper et al.^^ found that the 
critical current for given contacts does not depend on the 
orientation of the current flow above T 2.3 K. How- 
ever, below 2.3 K, a flnite difference continuously appears 
for supercurrents flowing in the positive or negative di- 
rection. This is in good qualitative agreement with our 
simple network model. 

6. Conclusions 

In this study, we investigated the phase diagram of the 
3K phase of eutectic Sr— 2Ru04-Ru samples, where su- 
perconductivity nucleates on the interface of Ru metal 
inclusions. For a single Ru inclusions, we find a sequence 
of two phases. At T* w 3 K, a time-reversal-symmetry- 
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Fig. 6. DiflFerence between the critical currents for opposite direc- 
tions of I: /(f and for critical currents in positive and negative 
directions, respectively. The plots have been taken using the pa- 
rameters, Iic/l2c = 4 and LIic/^o = 0.1 (solid line) and 0.05 
(dashed line). 



conserving phase (A-phase) nucleates on the interface. At 

a lower temperature T', a transition to a time-revcrsal- 
symmetry-breaking phase (B-phase) appears. This tran- 
sition is discontinuous (weakly first order) because the 
A- and B-phascs have different phase winding numbers 
at a closed (cylindrical) interface. In contrast to the A- 
phase, the B-phase carries spontaneous supercurrents at 
the interface due to the TRSB. These currents arc most 
likely responsible for the anomalous behavior of the crit- 
ical current, which is different for positive and negative 
current directions, as observed in an experiment below 
2.3 K.^^ Thus, we propose to identify the onset of this 
anomalous behavior as the transition between the A- and 
B-phascs. Further experiment on the quasiparticlc tun- 
neling that finds the onset of zero-bias anomalies in the 
tunneling spectrum below 2.4 K also indicates a phase 
transition within the 3K phase. This behavior also 
fits well within our phase diagram and will be discussed 



in detail elsewhere. 
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